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Abstract 

A two-equation k-co turbulence model has been devel- 
oped and applied to a quasi-three-dimensional viscous 
analysis code for blade-to-blade flows in turbomachinery. 
The code includes the effects of rotation, radius change, 
and variable stream sheet thickness. The flow equations 
are given and the explicit Runge-Kutta solution scheme is 
described. The k-co model equations are also given and 
the upwind implicit approximate-factorization solution 
scheme is described. Three cases were calculated: transi- 
tional flow over a flat plate, a transonic compressor rotor, 
and a transonic turbine vane with heat transfer. Results 
were compared to theory, experimental data, and to results 
using the Baldwin-Lomax turbulence model. The two 
models compared reasonably well with the data and sur- 
prisingly well with each other. Although the k-co model 
behaves well numerically and simulates effects of transi- 
tion, freestream turbulence, and wall roughness, it was not 
decisively better than the Baldwin-Lomax model for the 
cases considered here. 

Introduction 

A large percentage of computational fluid dynamics 
(CFD) analysis codes for turbomachinery use the Bald- 
win-Lomax turbulence model [1]. This was evident in the 
results of the blind test case for turbomachinery codes 
sponsored by ASME/IGTI at the 39th International Gas 
Turbine Conference held in The Hague in June of 1994. 
The results have not yet been published. Of the 12 partic- 
ipants, nine used the Baldwin-Lomax turbulence model, 
one used an algebraic mixing length model, and two used 
k-E models. One of the objectives of that test case was to 
investigate the effects of turbulence models. However, 
because of differences in grids, large variations between 
the computed solutions, and lack of experimental mea- 
surements in the boundary layers, it was not possible to 
draw any conclusions regarding the effect of turbulence 
models. 
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The Baldwin-Lomax model is popular because it is 
easy to implement (at least in 2-D) and works fairly well 
for predicting overall turbomachinery performance. How- 
ever, the model has both numerical and physical problems. 
Numerical problems include awkward implementation in 
3-D, difficulty in finding the length scale [2], and slow 
convergence if the length scale jumps between grid points. 
Physical problems include a crude transition model and 
the neglect of freestream turbulence, surface roughness, 
and mass injection effects which are often important in 
turbines. These effects are sometimes added to the Bald- 
win-Lomax model using techniques developed for bound- 
ary layer codes [3, 4]. Physical problems also include 
poor prediction of separation [5], which is important in 
compressors, and underprediction of wake spreading [2]. 

A few researchers have used other turbulence models 
for turbomachinery problems. Choi et. al. have used the q- 
co model [6], Hah (who participated in the blind test case) 
used a k-e model [7], and Kunz and Lakshminarayana 
used an algebraic Reynolds stress k-e model [8]. Unfortu- 
nately none of these researchers have used a Baldwin- 
Lomax model in the same code for comparison. Ameri 
and Amone have compared the q-co, k-e, and Baldwin- 
Lomax models for turbine heat transfer problems [9, 10], 

Two papers have compared the k-co and Baldwin- 
Lomax models for turbomachinery problems, Bassi, et. al. 
examined a film-cooled turbine cascade [11], and Liu et. 
al. examined a low pressure turbine cascade [12]. Both 
papers compared the computed results primarily with 
experimental pressure distributions. 

In the present work the k-co model developed by Wil- 
cox [13] was incorporated in the author’s quasi- three- 
dimensional (quasi-3-D) turbomachinery analysis code 
[14]. The code includes the effects of rotation, radius 
change, and stream surface thickness variation, and also 
includes the Baldwin-Lomax turbulence model. The k-co 
model was chosen for several reasons. First, the effects of 
freestream turbulence, surface roughness, and mass injec- 
tion are easily included in the model [13]. Second, transi- 
tion can be calculated using the low-Reynolds-number 
version of the model [15]. Third, Menter has shown that 
the k-co model does well for flows with adverse pressure 
gradients [5 ,16]. Finally, the k-co model should behave 
well numerically since it avoids the use of the distance to 
the wall and complicated damping functions. 
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Figure 1. Quasi-three-dimensional stream surface for a 
compressor rotor. 

This paper describes the quasi-3-D flow equations and 
the explicit Runge-Kutta scheme used to solve them. The 
paper also gives the k-o equations written in a quasi-3-D 
form, gives the boundary conditions, and describes the 
implicit upwind ADI scheme used to solve the turbulence 
model equations. The model was tested on three cases and 
compared to the Baldwin-Lomax model and to experimen- 
tal data. The cases included a flat plate boundary layer 
with transition, a transonic compressor rotor, and a tran- 
sonic turbine vane. 


e = p[C v r+i(u 2 + v 2 )J 

is the total energy per unit volume, 

P = (y- 1) [e-f w 2 + v 2 )J 
is the pressure, 

jn _\dr_ , jm _ 1 dh_ 

r 7 dm ' h hdm 

are derivatives of the stream tube geometry, and 
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Quasi-3-D Navier-Stokes Equations 

The Navier-Stokes equations have been developed in 
an (m, 0) coordinate system as shown in figure 1. Here 
m is the arc length along the surface, 

dm 2 = dz + dr 2 (1) 

and the 0 -coordinate is fixed to the blade row and rotates 
with angular velocity 0 . 

The radius r and the thickness h of the stream surface 
are assumed to be known functions of m. The equations 
have been mapped to a body-fitted coordinate system, sim- 
plified using the thin layer approximation, and nondimen- 
sionalized by arbitrary reference quantities p 0 , c Q , and 

. The Reynolds number Re and Prandtl number Pr are 

defined in terms of these reference quantities. The final 
equations are given in [14] and are summarized below. 

fy + ^E + ajF-Re" 1 ^ = K ( 2 ) 

where 

<7 = J- X \p,pu, pvr, <?] 7 
K = r l [o, K 2 , 0, of (3) 

S = J~'[o,S 2 , S 3 , S 4 f 


K 2 = (pv 2 + />-Re 1 CT 22 )y + (p-Re *033 
The viscous fluxes are given by 
S 2 = +T le a i 2 

-S 3 = (t\ m G l 2 + lc \e<*22) r 


)t (9) 
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a z sb y/?/p is the speed of sound squared. Using Stokes’ 
2 

hypothesis, X = , the shear stress terms are given by 
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U and V are relative contravariant velocities 
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u = % m u + \ Q w 

v = + Tle^ 


( 12 ) 


where w is the relative tangential velocity, w = v - r© . 

The 0-metrics are scaled by 1/r and the Jacobian is 
scaled by rh . The metric terms are found using central 
differences and 




= J 


0 n ~ m r, /r 

r 0 ^ m \ /r \ 


(13) 


J= [rhim^Q^-m^Q^)] 
The effective viscosity is 

\x = H L + li T 


-l 


(14) 


where the molecular (laminar) viscosity |i L is evaluated 
using a power law function of the temperature, and the tur- 
bulent viscosity |i T is evaluated using either the Baldwin- 

Lomax model [1] or Wilcox’s k-co model [13, 15]. Minor 
modifications to the coefficients and blending functions 
used in the Baldwin-Lomax model are described in [2]. 


Boundary Conditions 

At the inlet the total pressure, total temperature, and 
tangential velocity component are specified and the 
upstream-running Riemann invariant based on the axial 
velocity is extrapolated from the interior. At the exit, three 
of the four conserved variables are extrapolated and the 0- 
averaged pressure is specified using the method described 
by Giles in ref [17]. Periodic boundaries between the 
blades are solved like interior points using a dummy grid 
line outside the domain. 

Multistage Runge-Kutta Scheme 

The flow equations are discretized using finite differ- 
ences and solved using an explicit Runge-Kutta scheme. 
A spatially-varying time step and implicit residual 
smoothing are used to enhance convergence. Details of 
the solution scheme used here are given in (18) and are 
described briefly below. 

The discrete equations are solved using the explicit 
multistage Runge-Kutta scheme developed by Jameson, 
Schmidt, and Turkel [19]. A four-stage scheme is used. 
For efficiency, physical and artificial dissipation terms are 
computed only at the first stage. The Baldwin-Lomax 
model is updated every five time steps. The k-co model is 
usually updated every two time steps with twice the At of 
the flow solution. 


The spatially-varying time step is calculated as the 
harmonic mean of inviscid and viscous components in 
each grid directional. 

Artificial dissipation consisting of blended second and 
fourth differences is added to prevent point decoupling 
and to enhance stability. Eigenvalue scaling, as introduced 
by Martinelli and Jameson [20] but modified by Kunz and 
Lakshminarayana [8], is used to weight the artificial dissi- 
pation in each direction. The scaling is based on a blend 
of the one-dimensional time step limits at each point. The 
artificial dissipation is also reduced linearly by grid index 
near the wall and wake centerline to minimize the effects 
on the boundary layer. 

The explicit four-stage Runge-Kutta scheme has a 
Courant stability limit of about 2.8. Implicit residual 
smoothing introduced by Jameson and Baker in [21] can 
be used to increase the time step, and hence the conver- 
gence rate, by a factor of two to three. On high aspect 
ratio grids the stability limit is dominated by the grid spac- 
ing in the finest direction, and it is sufficient to use implicit 
smoothing in that direction only. The stability analysis 
given in [21] is used to calculate the smoothing parameter 
required at each point, then the same Eigenvalue scaling 
used for the artificial dissipation is used to reduce or elim- 
inate the smoothing parameter in grid directions where it 
is not needed. The use of Eigenvalue scaling for both the 
artificial dissipation and implicit smoothing greatly 
increases the robustness of the numerical scheme. 


k-co Turbulence Model 

The k-co turbulence model was first postulated by 
Kolmogorov in 1942 and later independently by Saffman 
in 1970 (see Wilcox’s book [13] for references.) It has 
been under development by Wilcox for many years and is 
described in detail in [13]. The model solves two turbu- 
lence transport equations for the turbulent kinetic energy k 
and the specific dissipation rate co. The model has a basic 
formulation for fully turbulent flows that satisfies the law 
of the wall without knowledge of the distance to the wall 
or complicated near-wall damping terms. There is also a 
low-Reynolds-number formulation used for modeling 
transition [15]. Boundary conditions can be specified to 
simulate mass injection or surface roughness. 

Most of Wilcox’s development of the model used 
boundary layer codes, but recently Menter has shown sev- 
eral applications to Navier-Stokes codes [16]. Menter 
found that the model exhibited strong dependence on 
freestream values of co and proposed a somewhat ad hoc 
fix. In this work many of Menter’s suggestions for numer- 
ical implementation of the model have been used, but his 
fix for the problem of freestream dependence has not. 

The quasi-3-D form of the k-equation has been 
derived by writing the m- and 0-momentum equations in 
non-conservative form, multiplying each by its fluctuating 
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velocity component, and Favre averaging- The usual tur- 
bulence modeling approximations are made, i. e., the 
Boussinesq model is used for the Reynolds stress terms, 
pressure work, diffusion, and dilatation are all neglected, 
and turbulent dissipation is taken to be proportional to 
k x co . The production term is written in terms of the vor- 
ticity magnitude using Menter’s suggestion [16], Source 
terms that arise from the quasi-3-D equations are 
neglected. The co-equation is derived from the k-equation 
by dimensional considerations. Wilcox’s constants are 
used without modification. The final form of the model 
equations is as follows: 

df+Ud^q+Vd^-Re^-G = l -(P-D) (15) 

where 


p*= (9/100) F p , a= (5/9) (F a /F^) , and <x*= 
where 



5/18+ (Re T /Rp ) 4 
1 + (Re T /R p ) 4 

«0 + Re T /R co 

1 + Re T /R w 

<Xq + Re T /R t 
1 + Re-p/R^ 

1L 

\l L (0 


with a 0 = 1 / 10 , (Xq = p/3 = 1/40, 


( 22 ) 



(16) andRp = 8 , R w = 27/10, R* = 6 . 


4 t = a* 


P* 

co 


(17) 


The molecular plus turbulent diffusion terms G are 
written using the thin-layer approximation giving 


Boundary Conditions 

At the inlet the turbulence intensity Tu and turbulent 
viscosity (i T are specified. Then k and co are found from 


f< + T le) 

(V + CT'H-p^* 

J 

(p. + crp T ) 9^(0 


( 18 ) 


Menter’s form of the production terms is used [16]. 


P 

P 


where 


R e-1 ^« 2 -|*(a m « + ;a e v) 
a[oc*Q 2 -|a)(a OT « + 


(19) 


Q. = d v - - 3 e u + v— ( 20 ) 

m r 6 t 

is the vorticity. The destruction terms are given by 


D _ p’cojfc 
P [pco 2 _ 

The baseline k-co model has five coefficients: 


( 21 ) 


P = 3/40, p*= 9/100, o= 1 / 2 , c*= 1/2 , a = 5/9, 
and the trivial constant a* = 1 . 

The low-Reynolds-number model replaces three of 
the constants with the following bilinear functions of the 
turbulence Reynolds number Re T : 


k = 2 TU %2 

to = a*^ 

Pt 


(23) 


where a* = 1 for the baseline model or a* = for the 
low-Reynolds-number model. Substituting equation (22) 
for a* into equation (17) for p. 7 . gives a quadratic for 

Re T . The solution is 




(24) 

and co may be found from 



(25) 

P L Re T 

A turbulent length scale can be defined using (17) as 


U T = = pfk(^£) = p Jk l (26) 

The effects of varying the inlet values of co or / is dis- 
cussed with the results. 

On solid walls k-0, and co is set using Wilcox’s 
roughness model. 
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wall 
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m 
100 
k R 


k+<25 
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(28) 


Udy-U+^j-q^j) -U (q u -q ul j> 

Ifi- = i(LT±l«7I) 

The diffusive terms (18) can be written as 



(V + crV T )d n * 


8 k d ri k 

J 

(p + aii-j-^co 


g 


(33) 


and is the equivalent sand grain roughness height in 

turbulent wall units. For all results shown here was set 
to 5, giving a hydraulically smooth surface. 

To avoid numerical difficulties near leading edges 
where \i x becomes large, an upper limit was imposed on 

co using a boundary condition suggested by Menter [16]. 


10 6v 

CO = — X 

max Re pA;y2 


(29) 


where Ay is the grid spacing at the wall. 

k and co were extrapolated at the exit and treated as 
periodic across trailing edge wake cut lines and between 
blade rows. 


ADI Solution Scheme 

An alternating direction implicit (ADI) scheme was 
used to solve the k-co equations. An implicit scheme was 
chosen so that the equations could be updated less often 
than the flow equations without stability problems, and 
also because the k-co equations are dominated by complex 
source terms which are evaluated only once in an ADI 
scheme, but potentially at every stage of a Runge-Kutta 
scheme. 

Equation (15) may be written as 

d,q = -[t/a^ + Va^-Re-'^G-^/’-D)] (30) 

Using first-order backward Euler time differencing 
and linearizing the right hand side about the previous time 
step gives 

{/ + At [t/a^ + Vd^ - Re' 1 j^G' - i ( P ' -£>')]} A q 

= -At\ud % q + Vd^q-Rt-^G-UP-D)] 

V (31) 

where primes indicate Jacobians. The advective terms are 
approximated using first-order upwind differences. 


Since a* = a = 1/2, g k = = g , and the Jacobian 

G v ' is simply 


G'Aq 


ac 

dq 


A<7 = g 


d^Ak 


(34) 


The diffusive terms are approximated by second- 
order central differences. 


W = \ 1 (8 u + 1 + g iJ {k iJ * 1 - hp (35) 

+ (gij.i + gij) (*i,;-i- 

The difference approximations (32) and (35) are diago- 
nally dominant and have zero row sum, which according 
to Baldwin and Barth [22], makes the implicit operator an 
M-type matrix with a non-negative inverse. 

Menter’s linearization is used for the source terms 
[16]. The production terms are treated explicitly, i.e., 
P' = 0. The destruction terms are linearized using 


= P’o) p*/t 
P [ 0 2(3co 


(36) 


The term in the upper right comer p*& is the only cou- 
pling between the k and to equations. Here the term has 
been neglected and the equations are solved uncoupled 
from each other and from the flow equations. 

Equation (31) is solved using an approximate factor- 
ization. 

[' * 4 ( w « * 7)] [' * 

= -At\ud^q + Vd n q-R^ J -G v -hP-D)] 

P H (37) 

The destruction terms D' are included with the 
streamwise ^ operator. Treating the destruction terms 
implicitly improves the diagonal dominance but gives the 
implicit operator a nonzero row sum. 
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Figure 2. Theoretical and computed skin friction for a 
flat plate boundary layer. 


The Baldwin-Lomax and k-co models have both been 
coded fairly efficiently. A turbulent flow solution updat- 
ing the Baldwin-Lomax every time step model takes 1.53 
times as long as a laminar solution. Past experience has 
shown that it is sufficient to update the Baldwin-Lomax 
model every five time steps, reducing the CPU time to 
about 1.1 times that of a laminar solution. A solution 
updating the k-to equations takes 1.6 times as long as a 
laminar solution. Some cases have been successfully 
computed updating the k-co equations every five time steps 
(with five times the At of the flow solution), making the 
k-co model nearly as fast as the Baldwin-Lomax model. 
Other cases failed to converge unless the k-co equations 
were updated every other time step, and that strategy has 
been used for all results reported here. The net result is 
that a flow solution with the k-co model takes about 1.18 
times as long as a solution with the Baldwin-Lomax 
model. 

Results 


Flat Plate Boundary Layer 

A turbulent flat plate boundary layer test case was 
used to verify the code. The flow had a Mach number of 

0.3, a Reynolds number based on plate length of 10 7 , and 
a freestream turbulence of one percent. The grid had 127 
points stream wise with 12 points upstream of the leading 
edge, and was stretched exponentially to resolve the tran- 
sition region. There were 41 points above the plate with a 

wall spacing of y + = 2 to 4 . 

Figure 2 compares the computed skin friction using 
both the Baldwin-Lomax and k-co models to a 1/5 power 



Figure 3. Theoretical and computed velocity profiles for 
a flat plate boundary layer. 


law correlation for laminar flow and a 1/7 power law cor- 
relation for turbulent flow. The k-co results agree closely 
with the turbulent correlation and the Baldwin-Lomax 
results are slightly low. 

The Baldwin-Lomax model was run with the transi- 
tion model given in their original paper [1], that is, the tur- 
bulent viscosity was calculated for each grid line but set to 
zero for the whole line if the maximum was less than 
cmutm = 14 times the laminar viscosity. The Baldwin- 

Lomax model transitions at Re x = 3xl0 5 , while the k-co 

model transitions at Re x = 5xl0 5 . Both results are 
within the experimentally observed range. The k-co 
model transition point was found to be strongly dependent 
on the freestream turbulence level, but fairly independent 
of the freestream value of co. However, the transition point 
calculate for a turbine vane shown later was found to be 
very sensitive to the freestream value of co. The Baldwin- 
Lomax transition point was strongly dependent on the 
value of the cutoff parameter cmutm , suggesting that this 
parameter could be calibrated to simulate freestream tur- 
bulence effects. 

Figure 3 compares computed velocity profiles located 
at the end of the plate to Spalding’s composite sublayer/ 
law-of-the-wall profile. Both results agree closely, but the 
Baldwin-Lomax profile has a slightly stronger wake com- 
ponent. 

Transonic Compressor Rotor 

The transonic compressor rotor described by Suder et 
al. [3] was used to test the quasi-3-D effects in the model. 
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Figure 6. Computed and measured near wake Mach Figure 7. Computed and measured far wake Mach 
number profiles for a transonic compressor rotor. number profiles for a transonic compressor rotor. 


The rotor was tested experimentally at NASA Lewis 
Research Center using both laser anemometry and con- 
ventional aero probes. 

A section of the rotor at 70 percent span was ana- 
lyzed. The radius was specified as a line 70 percent of the 
way between the hub and shroud, and the stream surface 
thickness was specified as the local distance between the 
hub and shroud, normalized to one at the inlet. A meridi- 
onal view of the stream tube is shown in figure 1. A C-type 
grid was used with 319 points around the blade and 45 

points away from the blade. The grid spacing gave y + < 2 
over most of the blade. The calculations were run 2000 
iterations, which took about 3.25 minutes for the k-co 
model on the Cray C-90 called eagle at NASA Ames 
Research Center. 

Figure 4 compares contours of relative Mach number 
computed with the k-co model (the Baldwin-Lomax model 
gives identical contours) to contours measured experimen- 
tally using laser anemometry. The inlet Mach number is 
about 1.4. The flow passes through a weak upstream -run- 
ning wave system, then through a strong normal passage 
shock, and leaves the rotor at about Mach 0.8. The com- 
puted Mach number behind the normal shock is somewhat 
lower than the measured values. This is partly due to the 
assumed stream surface but may also be partly due to the 
inability of either turbulence model to capture the shock- 
boundary layer interaction correctly. 

One motivation for using a turbulent transport model 
was that algebraic models frequently fail to find the cor- 
rect length scale and thus give nonsmooth values of \x T 

[2]. Occasionally nonconvergence or instability problems 
can be traced to poor numerical behavior of the turbulence 


model. Figure 5 shows contours of | i T computed using 

the Baldwin-Lomax model (left) and the k-co model 
(right.) 

Figure 4 also shows that the computed results seri- 
ously underpredict the wake spreading. This is shown 
quantitatively in figures 6 and 7. Figure 6 compares com- 
puted and measured near-wake profiles about 0.28 chords 
downstream of the trailing edge while figure 7 compares 
far- wake profiles about 2 chords downstream. The two 
turbulence models give surprisingly similar results, espe- 
cially considering the erratic behavior of the Baldwin- 
Lomax model seen in figure 5. 

The computed wakes are both narrower and deeper 
than the measured wakes. Wilcox has shown in [13] that 
his model gives the best prediction for planar wake 
spreading as co^— »0. This limit corresponds to 

\i T — > o© , which seems unreasonable for an inlet value. 

Varying 0 )^ by five orders of magnitude had very little 

effect on the computed wake spreading. The ffeestream 
turbulence was set to three percent, and doubling it had lit- 
tle effect on the wake spreading. 

Transonic Turbine Cascade 

A transonic turbine vane tested by Arts, et al. [23] was 
computed as a third test case. The vane was tested experi- 
mentally in the Isentropic Light Piston Compression Tube 
Facility at the von Karman Institute. The facility has inde- 
pendent control over the exit Reynolds number Re 2 ^ s , the 

exit Mach number, M 2 ^ s , and the inlet turbulence inten- 
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Figure 8. Computed Mach number contours for the VKI 
turbine vane. 

sity Tu. Surface pressures were measured with static taps, 
and wake total pressure profiles were measured with a 
high-speed traversing probe. The vanes were initially at 
300 K and the freestream temperature was 415 K. 
Unsteady blade surface temperatures were measured dur- 
ing a run using platinum thin film gauges, then converted 
to heat fluxes using a one-dimensional semi-infinite-body 
model. 

For the computations a C-type grid was used with 383 
points around the vane and 49 points away from the vane. 

The grid spacing gave y + < 1.5 over most of the vane. 
Blade shapes and computed Mach contours for 
^2 is = 0-9 are shown in figure 8. The flow accelerates 
from a Mach number of about 0.15 at the inlet to about 0.9 
at the exit. 

Surface heat transfer was converged to plotting accu- 
racy in 3000 iterations in fully turbulent regions. Laminar 
parts of the flow took longer to converge, so all calcula- 
tions were run 5000 iterations. Solution times were about 
eight minutes per case on the Cray C-90 computer. A typi- 
cal residual history is shown in figure 9. The k-co calcula- 
tions converged monotonically, but with the Baldwin- 
Lomax model the maximum residual oscillated near the 
round trailing edge. 



Figure 9. Residual histories for the VKI turbine vane 
computations. 

Computed distributions of isentropic surface Mach 
number are compared to experimental data for M 2 is = 

0.875 and 1.02 in figure 10. The Baldwin-Lomax and k-co 
models give identical results and are not shown separately. 
The subsonic results agree very well with the experimental 
data. The transonic results slightly underpredict the Mach 
number on the rear (uncovered) part of the suction surface. 
All subsequent results are for the subsonic case. 

Computed wake profiles located 43 percent of axial 
chord downstream of the trailing edge are compared to the 
experimental data (digitized manually from [23]) in figure 
11. Again the computed wakes are narrower and deeper 
than the measured wakes, but here the k-co results are 
slightly better than the Baldwin-Lomax results. 

Figures 12 - 15 show the effects of various parameters 

on surface heat transfer coefficient H [w/(m 2 K )]. In 
each figure the abscissa is the arc length S [mm] along the 
vane surface. Figure 12 compares computations using the 
baseline k-co model, the low-Reynolds-number k-co model, 
and the Baldwin Lomax model to experimental data for 

Re 2,/j? = lx 10 6 and Tu = 4 percent Triangles show the 

experimental data. The baseline k-co solution is fully tur- 
bulent on the suction surface giving high values of H. The 
pressure surface has a highly favorable pressure gradient 
and acts laminar over part of the chord, transitioning near 
S = -20. The low-Reynolds-number k-co solution remains 
laminar on the pressure surface and transitions near the 
measured transition point on the suction surface; however, 
the transition point was forced by choice of ja T ^ n , as dis- 
cussed later. The Baldwin-Lomax solution agrees closely 
with the low-Reynolds-number k-co solution. As dis- 
cussed with figure 14, the laminar parts of the flow have 
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Figure 10. Computed and measured distributions of 
isentropic Macn number for the VKI turbine vane. 



Fig ure 1 1 . Computed and measured total pressure profiles 
0.43 chords behind the VKI turbine vane. 



pressure surface S (mm) suction surface 

Figure 12. Surface heat transfer coefficient predicted by 
three turbulence models. 



Figure 13. Effects of inlet turbulent length scale or 
viscosity on predicted heat transfer coefficient 




pressure surface S (mm) suction surface 

Figure 14. Effects of inlet turbulence intensity on heat 
transfer coefficient 


pressure surface S (mm) suction surface 

Figure 15. Effects of Reynolds number on heat transfer 
coefficient 
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augmented heat transfer due to freestream turbulence that 
none of the models predict. 

Figure 13 shows the effects of | i Tin or the corre- 
sponding length scale l in on H for the same flow condi- 
tions. The inlet turbulent viscosity was varied by about 3 
orders of magnitude to produce solutions that ranged from 
fully laminar to almost fully turbulent. Only a small range 
of values of |i T ^ gave transition near the measured loca- 
tion. Corresponding turbulence length scales are shown 
on the figure. The length scale that gives the best transi- 
tion location, = 3.5x10 2 , is about 1/1600 times the 

pitch or about 7 times the grid spacing at the wall. All 
subsequent calculations were run with this length scale. 
Although this strong dependence on \ij^ is disconcerting 

since appropriate values are not known in advance, Wilcox 
points out that transition in real flows is not simply a func- 
tion of Tu, but is also frequency dependent [15]. He also 
suggests that two coefficients in the low-Reynolds-number 
model could be adjusted to better match other flows, 
although this has not been attempted here. 

The effects of freestream turbulence intensity Tu are 
shown in figure 14. The effect on suction surface transi- 
tion location is modeled reasonably well by the low-Rey- 
nolds-number k-co model. The experimental data shows a 
strong augmentation of heat transfer at the leading edge 
and on the pressure surface as Tu is increased. Although 
the low-Reynolds-number k-co model depends directly on 
Tu, in laminar regions the model gives values of turbulent 
viscosity that are much too small to affect the heat transfer. 
Boyle has added algebraic correlations to the Baldwin- 
Lomax model to simulate freestream turbulence effects 
[4]. These correlations could be added to the k-co code, 
but they do involve the distance from the wall. 

The effects of Reynolds number Re 2 ^ s are shown in 
figure 15 for Tu = 4 percent. The data shows a large 
increase in heat transfer with Re 2 ^ s . The k-co results show 

qualitative agreement, but the magnitude of the heat trans- 
fer is underpredicted due to the failure of the model to pre- 
dict freestream turbulence effects. Effects on suction 
surface transition location are overpredicted, and predicted 
transition is too abrupt. The data shows transition on the 
pressure surface at the highest Reynolds number that is not 
predicted, although a small change in does cause 

pressure surface transition, as shown in figure 13. 

Concluding Remarks 

Wilcox’s k-co turbulence model has been added to a 
quasi-3-D Navier-Stokes analysis code for turbomachin- 
ery. The code includes the effects of rotation, radius 
change, and stream sheet convergence, and also included 


the Baldwin-Lomax turbulence model. The quasi-3-D 
flow equations and boundary conditions were described. 
An explicit multistage Runge-Kutta scheme with spa- 
tially-varying time step and implicit residual smoothing 
was used to solve the flow equations. The quasi-3-D form 
of the k-co model equations and boundary conditions were 
also described. An upwind implicit ADI scheme was used 
to update the turbulence model equations uncoupled from 
the flow equations. The numerical scheme was quite 
robust, but about 18 percent slower than the Baldwin- 
Lomax model. 

Calculations were made for three test cases: a flat 
plate boundary layer with transition, a transonic compres- 
sor rotor with significant quasi-3-D effects, and a transonic 
turbine vane. The flat plate calculations agreed very well 
with theory for both turbulence models. Transition predic- 
tions were reasonable for both models and suggested that 
the Baldwin-Lomax transition model could be calibrated 
to simulate free stream turbulence effects. The compres- 
sor rotor calculations showed very close agreement 
between the two turbulence models, but both models 
failed to capture the measured wake spreading. The tur- 
bine calculations showed very good agreement with mea- 
sured surface pressures for both turbulence models. 
Predicted wake profiles were thinner and deeper than mea- 
sured profiles, although the k-co model gave marginally 
better results. The Baldwin-Lomax model, the baseline k- 
co model, and the low-Reynolds-number k-co model were 
compared for heat transfer calculations. The Baldwin- 
Lomax model did reasonably well considering the simple 
transition model used. The low-Reynolds-low-Reynolds- 
numnumber k-co model showed a high sensitivity to inlet 
values of co, expressed as an inlet turbulent viscosity or 
length scale which are not generally known. The k-co 
model was able to capture the effects of inlet turbulence 
intensity on transition but not on augmentation of heat 
transfer in laminar regions. It may be possible to model 
this effect with a simple algebraic model. Effects of Rey- 
nolds number were predicted qualitatively. 

The k-co model exhibited some attractive numerical 
properties, but for the cases considered here, predictions 
were not decisively better than those made with the Bald- 
win-Lomax model. Other test cases may identify areas 
where the k-co model is significantly better than the Bald- 
win-Lomax model. Future work will extend the k-co 
model to three dimensions where algebraic models are 
poorly defined and difficult to implement. 
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